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Abstract. We solve numerically the equations of motion for the col- 
lapse of a shell of baryonic matter, made of galaxies and substructure of 
1O 6 M0-=-1O 9 -Mq, taking into account the dynamical friction and the pa- 
rameters on which it depends: the peaks' height u c , the number of peaks 
inside a protocluster Ntot multiplied by the correlation function £(r) eval- 
uated at r = 0, the filtering radius Rf and the core radius of the cluster of 
galaxies, tq. We show how the collapse time of a shell of baryonic matter 
depends on those parameters. 



1. Introduction 

The problem of the formation and evolution of clusters of galaxies has been 
one of the crucial topic over the last years (see e.g. Colafrancesco et al. 1989, 
Kaiser 1993, Antonuccio-Delogu 1992, Colafrancesco et al. 1995 and Ryden & 
Gunn 1987). It is well known that the formation of cosmic structures is strictly 
related to the evolution of the density perturbations: in the present paradigma 
of structure formation, it is generally assumed that cosmic structures of size ~ 
R form preferentially around the local maxima of the primordial density field, 
once it is smoothed on the filtering scale Rf. These linear density fluctuations 
eventually evolve towards the nonlinear regime under the action of gravitational 
instability; they detach from the Hubble flow at turn around epoch t m , given 
by: 



tm 



32Gp b 



1/2 

(l+zf 2 (1) 



where p b is the mean background density, z is the redshift and 5 is the mean 
overdensity within the nonlinear region. After the turn around epoch, the fluc- 
tuations start to recollapse when their overdensity reaches the value 5 = 1. Since 
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the density field depends on the matter that dominates the universe, the mean 
characteristics of the cosmic structures depend on the assumed model. In this 
context the most succesful model is the biased Cold Dark Matter (hereafter 
CDM) (see e.g. Liddle & Lyth 1993) based on a scale invariant spectrum of 
density fluctuations growing under gravitational instability. In such scenario the 
formation of the structures occurs through a " bottom up" mechanism. 
A simple model that describes the collapse of a perturbation of density is that of 
Gunn & Gott (1972, hereafther GG72). This is in contradiction with the predic- 
tions of CDM models. It is well known that in a CDM Universe, an abundant 
production of substructures during the evolution of the fluctuations is predicted. 
The presence of substructure is very important for the dynamics of collapsing 
shells of baryonic matter made of galaxies and substructure of 1O 6 M0^-1O 9 M , 
falling into the central regions of a cluster of galaxies. In presence of substructure 
it is necessary to modify the equation of motion: 

fr _GM_ 

dt 2 r 2 (t) { ' 

since the graininess of mass distribution in the system induces dynamical fric- 
tion that introduces a frictional force term. In a material system, the gravi- 
tational field can be decomposed into an average field, Fo(r), generated from 
the smoothed out distribution of mass, and a stochastic component of the field, 
F stoc h(r), generated from the fluctuations in number of the neighbouring parti- 
cles. The stochastic component of the gravitational field is specified assigning 
a probability density, W(F), (Chandrasekhar & von Neumann 1942, hereafter 
CvN42). In an infinite homogeneous unclustered system IU(F) is given by the 
Holtsmark distribution (CvN42) whilst in inhomogeneous and clustered systems 
W(F) is given by Kandrup (1980) and Antonuccio-Delogu & Atrio-Barandela 
(1992, hereafter AA92) and respectively. The stochastic force, F stoc h, in a self- 
gravitating system modifies the motion of particles as if it was a frictional force. 
In fact, a particle moving faster than its neighbours produces a deflection of 
their orbits in such a way that the average density is greater in the direction 
opposite to that of motion, causing a slowing down in its motion. Adopting the 
notation of GG72 (see also their eqs. 6 and 8) and remembering that T c q/2 is 
the collapse time in the absence of dynamical friction (GG72), one can write: 

- 1/2 

where p& is the critical density at a time U and pi is the average density inside 
Ti at U. The equation of motion of a shell of baryonic matter in presence of 
dynamical friction, using the dimensionless time variable r = can be written 
in the form: 

d 2 a 4irGpcj(l + 6i) 2 da 

d^ = aJ(t) C ° " V { ) 

(Antonuccio-Delogu & Colafrancesco 1994, hereafter AC94) where 5% is the over- 
density within ri and n is the coefficient of dynamical friction. Supposing that 
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there are no correlations among random force and their derivatives, we have: 



_ Jd*FW(F)F*T(F) 

v ~ W) ( j 

(Kandrup 1980), where T{F) is the average "duration" of a random force im- 
pulse of magnitude F, W(F) is the probability distribution of stochastic force 
(which for a clustered system is given in eq. 37 of AA92). 



2. The collapse time 

In a previous paper (Del Popolo & Gambera 1996, hereafter DPG96), one of us 
showed how the expansion parameter a(r) depends on the dynamical friction, 
solving eq.(Q) numerically. DPG96 assigned some values to r] and then solved 
eq.(Q) through a numerical method not taking into account the parameters on 
which 7) depends. In this paper, we examine how the dynamical friction coeffi- 
cient r\ varies according to the parameters and how the collapse time depends on 
them. We solved eq.(|5|) numerically and the other equations on which it depends 
for a outskirts shell of baryonic matter with 5 = 0.01 inside the spherical regions 
(protocluster), for different values of v c , Rf, ro and 3, where: 

• v c is the peaks' height; 

• Rf is the filtering radius; 

• ro is the parameter of the power-law density profile; Theoretical work 
(Ryden 1988) suggests that the density profile inside a protogalactic dark 
matter halo, before relaxation and baryonic infall, can be approximated 
by a power-law: 

p(r) = *g (6) 

where p ~ 1.6 on a protogalactic scale. 

• 3 is the product iVf t ■£(()) where N to t is the total number of peaks inside a 
protocluster and £(0) is the correlation function calculated in r = 0. AA92 
have demonstrated that in the hypothesis m av 3> 1M , where m av is the 
average mass of the subpeaks, the dependence of the dynamical friction 
coefficient on N to t and £(r) may be expressed as a dependence on a single 
parameter that we define as: 

3 = N tot ■ £(0) (7) 
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Figure 1. Collapse time r of a shell of matter made of galaxies 
and substructure when dynamical friction is taken into account, versus 
v c We assume a core radius of tq = \h~ l Mpc and a filtering radius 
Rf = 0.74/i -1 Mpc. Open circles: H = 10 3 ; filled circles: H = 10 2 . 



After having determinated rj solving numerically eq. (|5|), we get r as a 
function of v c , Rf, ro and S solving eq.(Q). We perform these calculations for 
different set of values of v c , Rf, ro and 3; the results are shown in figs. 1^4. 
Before commenting upon the figures, we want to remark that the dependence of 
r on 6 is qualitatively shown in fig. 5 of AC94. We observe that for 5 > 1CP 2 
the collapse time in presence of dynamical friction is always larger than in the 
imperturbated case but the magnitude of the deviation is negligible for larger 5, 
whilst for 6 < 10~ 2 the deviations increase steeply with lower 5. Then, having 
considered 5 = 0.01, the estimation we get for r in §4 must be considered as 
a lower limit. In figure 1 we show the collapse time in presence of dynamical 
friction, versus the peaks' height, for different values of S. In this picture, we 
show how r grows for larger values of v c and for larger values of S. Similarly, 
in fig. 2 we note how r increases for larger values of v c and of Rf. The slope 
of the curves confirm our prevision on the behaviour of the collapse of a shell 
of baryonic matter falling into the central regions of a cluster of galaxies in the 
presence of dynamical friction: the dynamical friction slows down the collapse 
(as DPG96 had already shown) and the effect, as we are showing in figs. 1 and 
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Figure 2. Collapse time r of a shell of matter made of galaxies and 
substructure when dynamical friction is taken into account, versus v c 
We assume a core radius of tq = lh~ 1 Mpc and a fixed correlation 
S = 10 3 . Open circles: Rj = d.l^h~ l Mpc; filled circles: Rj = 
0.65/i _1 Mpc; crosses: Rf = 0.55h _1 Mpc. 



2, increases as E, Rf, v c grow. Here we want to remind that we are consider- 
ing only the peaks of the local density field with central height v larger than a 
critical threshold v c . This latter quantity is chosen to be the threshold at which 
fpeak iy > Vc) « lav where r pea k is the typical size of the peaks and l av is the 
average peak separation (see also Bardeen et al. 1986). 

In figs. 3 and 4 we show how the collapse time varies with the core radius 
of the protocluster ro- Note how r grows as ro decreases: the smaller the core 
of the protocluster, the larger the time of collapse in the presence of dynamical 
friction; besides we show how this effect increases for larger values of both v c 
and E. 
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Figure 3. Collapse time r versus ro We assume a filtering radius 
Rf = 0.74/i _1 Mpc. and a total number of peaks of substructure S = 
10 3 . Open circles: u c = 3; crosses: v c = 2; filled circles: v c = 1.24. 



3. Conclusions and discussion 

In the first part of this work we have shown how the collapse time r of a shell 
of baryonic matter made of galaxies and substructure depends on some parame- 
ters. In figures 1 and 2 we can see how r grows when E or Rf or v c increases. It 
means that the effects of the presence of the dynamical friction should be more 
evident in the outer regions of rich clusters of galaxies. Besides, we show how 
the collapse time of an infalling shell increases with decreasing values of ro , and 
becomes very large for ro < 2Mpc (see fig. 3). Then the slowing down of the 
collapse of an outer shell, of baryonic matter, within a cluster of galaxies owing 
to the dynamical friction is more remarkable in the clusters with core of little 
dimension. 
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Figure 4. Collapse time r of a shell of matter made of galaxies and 
substructure when dynamical friction is taken into account, versus ro 
We assume a filtering radius Rf = 0.74/i _1 Mpc and a peaks' height 
v c = 3. Open circles: S = 10 3 ; filled circles: E = 10 2 . 
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